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ABSTRACT 

We extend and investigate the spheroidal model of interstellar dust grains used to 
simultaneously interpret the observed interstellar extinction and polarization curves. 
We compare our model with similar models recently suggested by other authors, study 
its properties and apply it to fit the normalized extinction A(\)/Ay and the polarizing 
efficiency P(X)/A(X) measured in the near IR to far UV region for several stars seen 
through one large cloud. We conclude that the model parameter ft being the angle 
between the line of sight and the magnetic field direction can be more or less reliably 
determined from comparison of the theory and observations. This opens a way to 
study the spatial structure of interstellar magnetic fields by using multi-wavelength 
photometric and polarimetric observations. 

Key words: ISM: clouds - polarization - dust, extinction. 



1 INTRODUCTION 

Modelling of the wavelength dependencies of interstellar ex- 
tinction A(X) and polarization P(X) along with analysis of 
infrared bands is the main sour ce of information about the 
properties of interste llar grains |Drainell2003l ; IWhittet|[20ul ; 
I Voshchinnikovll2004 ) ■ 

Extinction data were often used to constrain 
the size distributio n of carbonaceous and s ilicate 
b are grains (e.g. iMathis. Rumpl fc Nordsieckl 1 19771; 



grains with the shape of these axisymmetric particles be- 
ing characterized by the only parameter — the ratio of 
the major to minor semi-axis a/b. The opti cal prop erties 
of spheroids have long been of intere s t (e.g., [Martini Il978l ; 
Rogers fc Martini Il979l: lOnakal [l980l; IDraine fc Led 1 1984 



Zubko. Krelowski fe Wegnerl 1 19961 ; Weingartner fc Draind 
T hese data were also modeled with com posite 



particles |Vaidva et al.ll200ll ; IVoshchinnikov et~al 2006; Iati 
let al.ll2008l ). In all these studies the grains were assumed to 
be spherical, i.e. their shape was not involved. 

To get information about the grain shapes as well 
as the ambient magnetic fields one needs to consider in- 
terstellar polarization. Interpretation of its wavelength de- 
pendence includes calculations of the extinction cross sec- 
tions of rotating p artially aligned non-spher i cal particles. In 
early studies (e.g. iHong fc Greenbera Il98dl ; IVoshchinnikovl 
1 19891 ) an unphysical model of infinitely long cylinders 
was appliecQ. More realistic is the spheroidal model of 
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1 iLi fc Gre cnbcre ( 199^) also considered finite cylinders for mod- 
elling interstellar extinction and polarization. 



Voshchinnikov fc Farafonovll99a ; lKim fc Martini 19951 ). The 
growth of such g rains in the interstella r med ium (ISM) has 
been discussed bv lStark. Potts fc Diverl (2006 ) , their absorp- 
tion and s catter ing properties (in the Rayleigh limit) b y 
iMin et all (|2006h and mineralogy bv lLi. Zhao fc Lil (|2007l ). 

Over the past decade spheroids have found wide 
application s in st udies of opti c s of cometary dust - 
iDas fc Sen! d200rj). iMoreno et all \200H ), circu mstellar dust 
' Wolf et al l (120021). erpla netary dust — lLasue et al.l 



J2007I). iGreenberg fc Lil dl996l) used prolate sphero ids and 



iLee fc Dra inc ( 1985) Jrlildebrand fc Dragovanl (| 19951 ) oblate 
ones to interpr et the silicate and ice po l arizat ion features of 
the BN object. IDraine fc Allaf-Akbaril (|2006l ) predicted ha- 
los due to X-ray scattering by oblate particles, detection of 
which would serve as a test for dust models. 

The spheroidal model of cosmic dust grains is particu- 
larly promising for interpretation of the interstellar polariza- 
tion and extinction data. A deta iled review of early works 
on the subject can be found in IVoshchinniko-vj (|2004|). so 
we co nsider recent applications of the model. Gupta et al.l 
(2005) well reproduced the average Galactic extinction curve 
in the region 0.3 — 9.5 ^m" 1 with silicate and graphite 
oblate spheroidal (a/b = 1.33) particles rotating in a plane 
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and having the size distribution n(rv) oc 



with rv 



0.005 — 0.25 /im where rv is the radius of a volume equiva- 
lent sphere. They found that the shape of grains did not af- 
fect the abundance problem, but essentially changed the po- 
larization efficiency factors at different particle orientation 
angles. It should be noted that interpretation of the extinc- 
tion data alone is usually made within the sph erical model o f 
dust grains (see for details the rece nt review of Drainej|2009) . 
IDraine fe Allaf-Akbaril (120061 ) and iDraine fc Fraissd (120091 ) 
fitted both the average Galactic extinction and polarization 
curves in the range 0.4 — 9.5 /im" 1 applying picket-fence ori- 
ented silicate and graphite oblate spheroids with the aspect 
ratio in the interval a/b — 1 — 2 with specific size distribu- 
tions. The grain alignment degree was size dependent and 
showed a sharp increase at rv ~ 0.1 fim. Voshchinnikov 
& Das (2008; hereafter VD08) studied the wavelength de- 
pendence of the ratio of the linear polarization degree to 
extinction (so called polarizing efficiency) from the ultravio- 
let to near-infrared. They found that the wavelength depen- 
dence of P(X)/A(X) was mainly determined by the particle 
composition and size whereas the absolute values of this ra- 
tio depended on the particle shape, degree and direction of 
alignment. 

In this paper we investigate further advantages of the 
spheroidal approach to interstellar grains. In contrast to 
other works we calculate the optical properties of spheroids 
of all sizes exactly even at far UV wavelengths, consider 
ensembles of particles with various degree and direction of 
alignment, compare different approaches to extinction and 
polarization curve modelling, treat the observational data in 
a new way (the wavelength dependence of the polarizing ef- 
ficiency is fitted instead of the polarization curve), apply the 
theory to observational data for concrete stars, and discuss 
a possibility to determine the direction of interstellar mag- 
netic fields from the data. In Section 2 we describe details of 
our model and compare it with the models used earlier. In 
Section 3 the results obtained from our model are compared 
with observations of seven stars for which the line of sight 
intersects one main cloud, keeping in view that "mean ex- 
tinction curves as well as mean total-to-selective extinction 
ratios are possibly averaged over several i nterstellar cloud s 
of possibly different physical parameters" (|KreIowskj|i989h . 
Section 4 contains our conclusions. 



2 THEORY 

Our modelling is based on the simple standard assumptions 
concerning dust grain materials, size distribution, shape and 
alignment. Homogeneous spheroids considered are charac- 
terized by their type (prolate or oblate), the aspect ratio a/b 
and a size parameter. As the latter we choose the radius rv 
of a sphere whose volume is equal to that of a non-spherical 
particle, i.e. r v = ab 2 for prolate spheroids and r v = a 2 b 
for oblate ones. 



2.1 Dust grain materials 

Interstellar dust mainly consists of five heavy elements: C, O, 
Mg, Si an d Fe, which are locked in carbonaceous and silicate 
particles (| Jones!l2000l ; iDrainef 2009 ) . We consider sub- micron 



grains from amorphous carbon (AC1) and amorphous sili- 
cate (astronomical silicate, astrosil). In order to reproduce 
the 2 175 A absorption feature small graphite spheres with 
radius r gla — 0.02 fim are also involve ql. The optical con- 
stants of t he materials ar e take n from iRouleau fc Martin] 
(|l99lh and lLaor fc Drainel (|l993l ). 



2.2 Size distribution 

Interpretation of the interstellar extinction is usually aimed 
at reconstruction of the dust grain size distribution consis- 
tent with the cos mic abundances and s ome other constraints 
(see discussion in I Voshchinnik ov 2004). 

We choose a power-law size distribution 

n(r v ) oc r~ q , (1) 

which wa s derived by Mathis. Rurrrpl fc Nordsieckl (|l977T l 
(see also IDraine fc Lea 1 19841 ) from minimization of x 2 - 
statistic. The distribution has three parameters: the lower 
(7V, m in) and upper (rv,msx) limits and the power index 
q. It was widely used in modelling of interstellar extinc- 
tion and radiative transfer in dusty objects. The average 
Galactic extinction curve can be reproduced with the so- 
called standard MRN mixture — an ensemble of carbona- 
ceous (graphite) and silicate spheres (in nearly equal pro- 
portions) with the parameters: q — 3.5, rv,min ~ 0.005 )im 
and 7V,max ~ 0.25 fim. 



2.3 Cross section averaging 

Let us consider non-polarized light passing through a cloud 
that contains three populations of grains: rotating spheroids 
from amorphous silicate (Si), amorphous carbon (C) and 
graphite (gra) spheres. For a line of sight, extinction (in 
stellar magnitudes) and linear polarization (in percentage) 
produced by the cloud can be written as 

•A(A) = 1.086A r d(C cx t)A , P(A) = Af d (C po i)Al00%. (2) 

Here, Ad = Ac + Asi + A gra is the total dust grain column 
density, (C ex t)x and {C po i)\ are the extinction and polar- 
ization cross sections, respectively, averaged over the grain 
populations, i.e. 

(C'ext)A = ^CcCcxt,c(A)+/CsiC'cxt,Si(A)-|-/Cg ra Cext,gra(A) , (3) 

(C P oi)a = /CcC po i i c(A) + A^siC'poi,si(A) , (4) 

where C cx t, g r a (A) = 7r7- gra Q ox t, g r a (A) and Ki = Ni/N d (i 
( '. Si) is the relative column density of the i-th population. 
It is evident that 

K,c + fcs\ + /C gra = 1. 



sections over the size distribution and grain 1 
l|Hong fc Greenberdll980l : IVoshchinnikovill989l ) 



The values of C ex t,pol,j are obtained by averaging of the 
cross sections over the size distribution and grain orienta- 
tions _ 

<?«m(a) = (§) y y y y ^q^w 

x m(rv) /(£, /3) dipduidfidrv 



7r/2 fir/2 pir/2 



(5) 



2 Note that graphite is not the only material consider ed as the 
carrie r of the the 2175A absorption feature (see, e.g., iLi et all 
|200S|) 
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Figure 1. Geometrical configuration of a spinning and wobbling 
prolate spheroidal grain. The major (symmetry) axis of the par- 
ticle O1O2 is situated in the spinning plane NO1O2 which is 
perpendicular to the the angular momentum J. The direction of 
light propagation k is parallel to the Z-axis and makes the angle 
a with the particle symmetry axis. 



Cpol,i(A) — 



tt/2 



tt/2 



TV.i Qpol,i(A) 



x n,i(rv) /(£, ft) cos 2ij) dipdu)df3drv ■ 



(6) 



Here, ft is the precession-cone angle for the angular mo- 
mentum J which spins around the direction of the mag- 
netic field B, ip the spin angle, uj the precession angle 
(see Fig. [I). The quantities Q cxt = (Q™ + Q™)/2 and 
Q P oi = (Q™ — Q2t)/2 are the extinction and polariza- 
tion efficiency factors for non-pol arized incident radiation 
|Voshchinnikov fc Farafono"vlll993l ). 



. cos 2<3> 1 , 

cos 2-0 = \ L 1 " 008 f sm 

cos 2$ 



2 cos ip COS 

7— 



prolate 
oblate, 



and 
cos 6 

cos 2$ 



cos S7 cos /3 + sin SI sin /3 cos cj, 



2 sin 2 /3 sin 2 w 



sin # cos tp 
a/1 - sin 2 



prolate 
oblate. 



(7) 

(8) 
(9) 

(10) 



Here, S7 is the angle between the line of sight and the mag- 
netic field direction (0° ^ Q «S 90°). The value SI = 90° 
corresponds to the case when the particle rotation plane con- 
tains the light propagation vector k, which gives the max- 
imum degree of linear polarization. For S7 = 0°, the light 
falls perpendicular to the particle rotation plane and from 
symmetry reasons the net degree of polarization produced 
is zero. 

The efficiency factors were calculated by using a so- 
lution to the light scattering problem for a homogeneous 



spheroid obtained by the me thod of separation of vari- 
ables in spheroidal coordinates (jVoshchinnikov fc Farafonovl 
Il993h and a mo dern approach to computations of s pheroidal 
wave functions l|Voshchinnikov fc Farafonovl 1200^ ). In con- 
trast to all works where th e spheroidal model was recently 
appli ed (|Gupta et all 120051 ; iDrainel 120091 ; iDraine fc Fraissel 
120091 ) we computed the optical properties of the particles of 
all sizes exactly even in the far UV spectral region. 



2.4 Grain alignment 

Interstellar non-spherical grains should be partially aligned. 
The "picket fence" or perfect rotat ing alignment often 
used in modelling of po larization (e.g. iKim fc Martiniri995l ; 
IDraine fc~F raissc 2009) is a crude approximation giving a 
too high polarizing e fficiency in compa rison with the empir- 
ical upper limit (e.g.. lGreenbersdll978l ) 

PmaxMv £ 3%/mag. 

Alt hough the theor y of grain alignment is actively dis- 
cussed (|Lazarianl 120071 ) it still has little practical signifi- 
cance. Therefore, we consider the classical alignment mecha- 
nism— sc^alledJrn^CTf^t^ayis-Greenstein (IDG) orienta- 
tion l|Davis fc Greensteinlll95ll ). The IDG mechanism is de- 
scribed by the function /(£, f3) dependi ng on the alignment 
parameter £ and the prece ssion angle /3 jHong fc Greenbera 
1 1980) ; IVoshchinnikovlll989l ) 



£sin/3 



(£ 2 cos 2 P + sin 2 j3) 3 / 2 



(11) 



The parameter £ depends on the particle size ry, the imagi- 
nary part of the grain magnetic susceptibility x" — xuid/Ta , 
where u>d i s the angular velocity o f a grain and x = 
2.5 x 10~ 12 (D avis fc Greensteinlll95ll ). gas density hh, the 
strength of magnetic field B, and temperatures of dust Ta 
and gas T g . By introducing a parameter 



5 = 8.23 x 10 2 



xB 2 



n H T K 2 T d 



one can get 
£ 2 = r v + S (T d /T g 



(12) 



(13) 



rv + So 

We assume that carbon and silicate grains have similar 
alignment but different size distributions. Thus, our model 
has 10 main parameters: ?v >m j n , ry m „ and q for carbon and 
silicate grains, relative density for carbon (ICc) and silicate 
(/Csi) grains, degree (So) and direction (Q) of grain align- 
ment. 



2.5 Comparison with other models 

There were just a few approaches to simultaneous mod- 
elling of the wavelength dependencies of interstellar ex- 
tinction and linear polarization. In most papers where an 
imperfect alignment of grains was involved (along with a 
power-law size distribution) the authors derived that parti- 
cles with ry Si 0.1/im should be poorly (randomly) oriented 
while those with ry <; 0.1/xm should be nearly perfectl y 
aligned (e.g., IKim fc Martini [l995l ; IDraine fc Fraissel 120091 ). 
Some physical grounds for such alignment of interstellar 
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grains come, e.g., from lMathisl (|l986l ). lGoodman fe Whittetl 
l| 19951 ) . iMartU l|l995h . 

As the actual com position of carbonaceou s grains in the 
ISM is still unknown (|Draine fe Fr aissc 2009), in interstel- 
lar dust modelling one often considered large (sub-micron) 
graphite grains instead of amorphous carbon ones. As po- 
larization produced by such graphite spheroids may be to a 
certain extent peculiar, the carbonaceous grains were often 
assumed to be either spherical or r andomly aligned (e.g., 
lMathis|[i"986l ; iDraine fe FraisseH2009T ). 

To compare these approaches with our one, we have 
made some calculations for three additional models: 

model Al - imperfectly DG-aligned silicate and graphite 
spheroids; 

model A2 - imperfectly DG-aligned silicate spheroids 
and randomly oriented amorphous carbon spheroids; 

model A3 - silicate and amorphous carbon spheroids 
with randomly oriented small particle and perfectly DG- 
aligned large particle. 

The grain size distributions and the population of small 
graphite particles are not varied, i.e. they are the same as 
in our basic model described in Sects. I2.1H2.4I 

When considering graphite spheroids in the Al model, 
we use the standard 2/3- 1/3 approximation which gives rea- 
sonably accurate results (jDraine fe Malhotralll993l 1. The ori- 
entation function for amorphous carbon grains in the A2 
model is f(P) = sin/3. In the A3 model this function is size 
dependent 



sin p for r v < rv.cut 
8{P) for r v > rv.cut 



(14) 



where 8(P) is the Dirac delta function, ry ]CU t a cut-off pa- 
rameter. 

So, we assume that all particles whose volume is smaller 
than that of a sphere with the radius 7v, cu t are randomly 
oriented, while other particles have perfect rotation align- 
ment, i.e. their angular moment vectors are parallel to the 
magnetic field dire ction. In our calculation s we choose ry jC ut 
to be 0.11 ^m (cf. IDraine fe Fraissdl2009T ). 

Some results obtained with the models considered are 
presented in Fig. [2] the fitting procedure used is described 
in the next Section. We see that all models can give very 
similar extinction curves that well fit the observational data. 
Only the Al model with large graphite particles provides the 
curve that differs a bit from others. Certainly, the extinction 
data are not enough to constrain the models. 

The model dependencies of polarizing efficiency shown 
in the right panel of the figure vary in form. In particular, 
the model with randomly oriented carbonaceous particles 
produces too small polarization in the visual and IR parts of 
spectrum. Apparently, in this case a better fit can be reached 
if one takes a more complex size distribution function or 
(and) increases the alignment efficiency of silicate grains. 
However, this complicates the model, therefore, we prefer 
our model having a small number of the parameters and 
simple assumptions. 



3 COMPARISON WITH OBSERVATIONAL 
DATA 

So far, most efforts in modelling of the optics of interstel- 
lar dust were directed at fitting of the average Galactic 
extinction curve (often with some additional constraints) 
by making use of spherical particles. Interpretation of the 
curves observed towards individual objects was very sel- 
dom, and that of the extincti on and pola rization curves was 
done only for a few objects. IH'inl (1987) applied the model 
of infinitely long coated cylinders with the imperfect DG- 
alignment to the interpretation of the wavelength depen- 
dence of interstellar extinction, linear and circular polariza- 
tion for t he star HP 204827 in the visual and IR parts of 
spectrum. iKim fe Martini l|l994h used infinitely long cylin- 
ders with the perfect DG-alignment (fi = 90°) to inter- 
pret the data for HD 25443 and HD 197770. Infinitely long 
cylinders with a silicate core and org anic refractory mantle 
were chosen bv lLi fe Greenbe"rel l|l998l ) for the interpretation 
of the extinction and po l arizat ion curves for HD 210121. 
So, IVoshchinnikov fe Dad (|2008l ) were first who utilized the 
spheroidal model to consider both the extinction and polar- 
ization wavelength dependencies observed for two stars. 



3.1 Sample of objects 

For application of our model we chose stars seen through 
a single cloud with the extinction and polarization curves 
known in a wide wavelength region. Initially, preference 
was given to five stars (HD 24263, HD 62542, HD 99264, 
HP 147933 and HP 19 7770) with measured UV polarization 
l| Anderson et al.ll 19961 ). These stars are not peculiar and not 
very distant. Later, two stars (HD 147165 and HD 179406) 
for which the polarization data in the visible were obtained 
by us was added to our sample. 

In this paper we discuss extinction and polarization of 
six stars (HD 24263, HD 99264, HD 147165, HD 147933, 
HD 179406, and HD 197770). Two objects (HD 62542 and 
HD 147933) were considered in our previous paper (VD08). 
Almost all the stars are not very far away from the galac- 
tic plane (Table [1} and are in dusty environment. One 
can assume that these stars are seen through single dust 
cloud s as they have weak r otation angle of linear polariza- 
tion (lAnderson et alj 19961) and narrow symmetric sodium 
Di and D? lines i|Zubko, Krelowski fe Wegnen[l996l ). A few 
comments on the stars should be made. 



HD 24263 is a double star (the secondary is a B9 star) seen 
through a large cloud AG6 located between the Local Bubble 
and Orion-Eridanus Super-bubbl e at distance ~150 pc and 
observed in Nal D and Call lines (|Genova fe Beckmanl f2003; 
I Welsh et al]|2005h . One als o detects the molecule s CH, CH+ 
and CN in the line of site (|Penprase et al.lll990l ). 

HD 99264 is located in the region of the Scorpios- 
Centau rs OB association and the Chameleon-Musca dar k 
clouds (jPlatais et al.lll998l ; ICorradi. Franko fe KnudeH2004T l. 
Consideration of high-resolution Nal D profiles and 
uvbyP photometry for over 6 nearby B stars allows 
ICorradi. Franko fe Knudel (|2004T ) to conclude that interstel- 
lar gas and dust in this region are distributed in two ex- 
tended sheet-like structures. The nearby feature located at 
the distance < 60 pc provides the colour excess E(b — y) ~ 
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Figure 2. Normalized extinction (left panel) and polarizing efficiency (right panel) for our basic model and three additional models with 
graphite spheroids (Al); randomly oriented amorphous carbon spheroids (A2); small randomly oriented silicate and amorphous carbon 
spheroids (A3) (see the text for mode details). Oblate spheroids with a/b = 3 and Q = 33°, for other parameters see the last line in 
Table [3] Observational data for HD 147933. 



0.05 while the second feature at 120-150 pc has E(b — y) f» 
0.2. 

HD 147165 (a Sco) is a beta Cephei type pulsating vari- 
able and a spectroscopic bin ary with the secondary being a 
B1V star (jNorth et al.ll2007h . It is a member of the Upper 
Scorpius subgroup within the Sco OB2 association related 
with reflection nebul a in the Ophiuchus cloud complex re- 
gion llMamaiekll200i ). 

HD 147933 is a component of the well-studied binary star 
p Oph AB. High-resolution KI, Nal, Call observations 
show very comple x profiles of atomic lines in its spectrum 
|Snow et al.ll2008l ). The star is associated with a strong re- 
flection nebula and a significant amount of dust (and gas) 
is assumed to lie behind it. 

HD 17 9406 (20 Aql) is a v a riable star behind a t ranslucent 
cloud (|Hanson et al.l Il992l ; iDobashi et. al] 120051 ). Though 
the absorption and emission lines have at least 3 compo- 
nents, the dominant one is seen in most of the atomic and 
all molecular lines. 

HD 197770 is an evolved spectroscopic eclipsing binar y star 
with both components being B2III stars l|Clavtonl 19961 ). The 
object is close to Cyg OB7 and Cep OB2 and probably lies 
on the edge of a star formation region including molecular 
clouds L1036 and L1049. It is as sociated with non-s tellar 
IRAS sources (see for more details [Gordon et al.lll99"^ ). 

So, the clouds observed in the selected lines of sight look 
to be rather different and hence their dust grain ensemble 
parameters may also differ. 



3.2 Observational data 

Basic information on the selected stars is presented in Ta- 
ble [1] The distance D, colour excess E(B — V), ratio of 
total extinction to the selective one Ry and extinction 
Ay for all stars exc e pt H P 24263 and 197770 are from 
iFitzpatrick fc Massal (|2007l ). with the data for these two 



stars being from IValencic. Clayton fe Gordon! l|2004h . The 
MK spectral types wer e taken fro m the catalogue of stellar 
spectral classifications (Skiff 2009). 

The photometric data in the standard visual - near in- 
frared bands were found in the gene ral catalogue of photo- 
metric data |Mermilliod et al] 1 19971 ) and recent papers re- 
ferred in the Simbad database (see Table [5J). When the JHK 
data were unavailable, the JHKg values from the 2MASS 
catalogue (|Cutri et al. 2003) were used. When the RI data 
were absent, we app lied with caution the data from the 
USNO Bl catalogue (|Monet et alj|2003l ). To derive the in- 
terste llar extin c tion in the line of sight the intrinsic colours 
from IStraizvs! l| 19921 ) were utilized. The UV extinction 
param eters were taken from Valencia Clayton fc Gordon! 
(|2004h and IFitzpatrick fc Massal (|2007l ) 

Uncertainties of UV extinction are mainly a result 
of errors in Ry and Ay values. Note that typical rel- 
ative errors of these quantites are 10-15% (see, e.g., 
IValencic. Clayton fc Gordon! 1 20041 ). An additional source of 
relatively large errors of extinction values in the visual and 
near-IR regions is uncertainty of the intrinsic colors (cf. , e.g., 
iBessell fc Brett||l988l ; rStraizvdli1^ : lDucati et al.ll200ll ). So, 
relative errors of extinction for individual objects are rather 
large (see, e.g., Fig. 2). 

Table [2] also contains references to the papers with 
polarimetric data. Additional polarimetric observations of 
HD 147165 and HD 179406 were performed by one of the 
authors (HKD) in 2007 at a 2m telescope of the Girawali 
Observatory (Pune, India). The IUCAA Faint Object Spec- 
trograph and Camera (IFOSC) was used in the polarimet- 
ric mode. The polarization parameters were measured in 10 
bands in the wavelength region about 0.37 - 0.60 /im and a 
standard package was applied for data analysis. 

Note that for two stars (HD 62542 and HD 197770), 
far-UV extincti on curves were de rived from IUE and FUSE 
observations bv lSofia et all |2005l ). For two stars (HD 147933 
and HD 197770), the important polarization data in the 
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Table 1. Target stars. 



Star 


I 


b 


D (pc) 


Sp 


Hy 




HD 24263 


182.1 


-34.9 


171 


B5V 


3.44 


0.72 


HD 62542 


255.9 


-9.2 


396 


B5V 


2.74 


0.99 


HD 99264 


296.3 


-10.5 


266 


B2.5V 


3.15 


0.85 


HD 147165 


351.3 


+17.0 


137 


B1III 


3.60 


1.48 


HD 147933 


353.7 


+17.7 


118 


B2.5V 


4.41 


2.07 


HD 179406 


28.2 


-8.3 


227 


B3V 


2.88 


0.89 


HD 197770 


93.9 


+9.0 


9 13 


B2V 


2.77 


1.61 



Table 2. Photometric and polarimetric data sources. 



Star 


Photometry^ 


Polarimetry^ 


HD 24263 


V04 (UV); M97 (UBV); M03 (RI); L05 (JHKL) 


A96 (UV & vis-nearlR) 


HD 62542 


S05 (farUV); F07 (UV); M97 (BV); M03 (RI); C03 (JHK S ) 


A96 (UV & vis-nearlR) 


HD 99264 


F07 (UV); M97 (UBVRI); C03 (JHK S ) 


A96 (UV & vis); S69 (UBV) 


HD 147165 


V04 (UV); W02 (UBVRI JHKLM) 


this work (vis); C66 (vis— nearlR) 


HD 147933 


F07 (UV); M97 (UBVRI); C03 (JHK S ) 


A96 (UV-nearlR) 


HD 179406 


V04 (UV); W02 (UBV JHKL); M03 (RI) 


this work (vis) 


HD 197770 


S05 (farUV); V04 (UV); M97 (UBV); M03 (RI); C03 (JHK S ) 


A96 (UV-nearlR) 



a M9 7 is for iMermilliod et alj il997h . W02 f or IWegnerl <2002l). C03 forlCutri et alj [|2003l1 , M 03 for iMonet et alj 
j2003ri , V04 for IValencic. Clayton fc Gordon] j2004h , L05 for lLarson fc WhittetJ 11200511 . S05 for lSofia et alj j2005ft . 



(ZUUdJ, V U4 tor valcncic, Clayton & 
F07 for lFitzpatrick fc Massal ^20071 ). 



b C66 is for lCovne fc Gehrelj < 1966ft . S69 for lSerkowski fc Robertson! j 196911 . A96 for I Anderson et alj dl996Tl . 



J, H and K bands are also available (|Wilking et al.lll980l; 
IWilking. Lebofskv fc Riekll 19821 ). 



3.3 Fitting procedure and parameter variations 

A general problem of dust optics modelling is in obtaining 
the absolute values of t he model parame ters. As it has been 
noted many times (e.g. . iGreenberdl 1978h . similar extinction 
occurs when a product of the typical particle size (r) and 
the particle refractive index does not change, i.e. 

(r) \m — 1| w const. 

The average Galactic extinction curve was successfully 
reproduced using quite different mixtures of particles 
with exponential dHong fc Greenberd Il980t ). power-law 
(|Mathis, Rumpl fc NordsieckHl977l)~ and more complicated 
grain size distributions fMathis 1996; Wcing artner fc Drainei 
l200ll l. Effects of variations of the distribution parameters 
(e.g., ?v,mm, rv.max, q for a power-law distribution) on 
extinction are well studied for spheres (see, for example, 
IVoshchinnikov fc Il'iij|l993h . 

Joint modelling of interstellar extinction and polariza- 
tion adds at least one parameter — particle shape — to the 
model. In the case of spheroids the particle shape is defined 
by the spheroid's aspect ratio a/b. Besides this, we can con- 
sider prolate ("cigars") and oblate ("pancakes") spheroids. 
If such particles are imperfectly oriented, the alignment pa- 
rameters (in our case, So and fi) very weakly affect extinc- 
tion. Note that the normalized extinction A\/Av essentially 
changes in the UV region when a/b varies (see Fig. 3 in 
VD08). 



In order to analyze importance of the model param- 
eters we reproduced the observational data available for 
HD 147933 with different particles. Some results are shown 
in Figs. [5] and |3 (see also Fig. 8 in VD08), the parameters 
of the fits are given in Table O Note that in VD08 the fit- 
ting was simply "chi-by-eye" while in this paper we apply 
minimization of x 2 . 

The fitting approach used included two steps; i) finding 
of the size distribution parameters and the relative column 
densities /Q from approximation of the extinction data using 
particles of given shape, and ii) a more accurate determina- 
tion of the alignment parameters So and Q from fitting of the 
wavelength dependence of polarizing efficiency P(X)/A(X). 

It is interesting that by varying the parameters rv, m i n , 
>V,max and q and the relative density of carbon (/Cc) and sil- 
icate (Ksi) grains one can reproduce the observational data 
with prolate as well as oblate particles of different aspect ra- 
tios. Note that after fixing the particle type (prolate/oblate) 
and a/b one can estimate only the values of /Cc, A^si, rv.min 
and q as the extinction and polarization curves is little af- 
fected by variation of rv,max- 

The alignment degree for three models presented in Ta- 
ble |3] is found to be slightly larger than the stan dard inter- 
stellar value (So ~ 0.2 nm. IVoshchinnikov| [l989). This can 
be a result of enhanced magnetic fields in the interstellar 
cloud near p Oph with which HD 147933 is related. The 
most exciting conclusion of our modelling is a close coinci- 
dence of the magnetic field directions — the values of Q for 
various fits differ less than 10°. We suggest that this model 
parameter can be more or less reliably determined from such 
modelling. 

Possibly, a better coincidence of the theory with the ob- 
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Table 3. Modelling results for HD 147933. 





amorphous 


carbon 






astrosil 




graphite 


shape; 






x 2 




ICc 






g 


/C si 


r V,min 




'/ 




a/b 


<5 (/im) 


n (dc g ) 


A/Ay 


P/A 


0.24 


0.03 


0.15 


1.5 


0.73 


0.08 


0.25 


1.5 


0.03 


pro; 5 


0.5 


35 


1.04 


24.6 


0.54 


0.10 


0.20 


2.0 


0.40 


0.01 


0.25 


1.3 


0.06 


pro; 3 


0.3 


12 


0.19 


8.2 


0.50 


0.10 


0.25 


2.0 


0.47 


0.07 


0.20 


2.0 


0.03 


obi; 3 


0.3 


33 


0.19 


2.2 




02468 02468 

A -1 , /xm _1 A -1 , /xm" 1 

Figure 3. Comparison of the normalized extinction (left panel) and the polarizing efficiency (right panel) observed for HD 147933 with 
results of modelling. The solid curves show the case of prolate (a/b = 3) spheroids having a power-law size distribution. Other model 
parameters are given in Table [3] 



servations could be obtained if we consider different align- 
ment parameters for silicate and the carbonaceous grains. 
Certainly, consideration of non-aligned carbon a ceous grains 
cann ot be ruled out ( iLi fc Greenberd |2002| ; IChiar et alj 
120051 ; see, however, Fig. [2}. 



3.4 Fitting results and discussion 

The normalized extinction A(X)/Ay and the polarizing effi- 
ciency P(X)/ A(X) were calculated with the theory of Sect. [5] 
and compared with the observational data available for the 
stars mentioned in Sect. 13.11 The results of this compari- 
son are plotted in Figs. EHU (see also Fig. 9 in VD08 for 
HD 62542) with the parameters of the fits being presented 
in Table [4] We start with comments on the individual ob- 
jects and then make general conclusions. 
HD 62542. The UV extinction curve has a very broad and 
weak 2175 A bump and a steep far-UV rise IjSnow fc McCalll 
2006). Our model suggests that the major part of grains are 
small carbon particles with a very narrow size distribution. 
HD 99264- Th e extinction curve for it is characterized by 
Rv of about 3 (|Fitzpatrick fc MassajbOOTl ). the UV extinc- 
tion curve is a bit "shallow" but not anomalous |Carnochanl 
1986). The polarization in an UV region and the UBV bands 
is well described by Serkowski curve with A max = 0. 55/jm 
l| Anderson et all 1 19961 ; [Martin. Clayton fc Wolj|l999l) . So, 



we have the case close to the average Galactic curves, but 
certainly there are not enough polarization data for sound 
comparison. 

HD 147165 and HD 179406. Extinction in these lines of 
sight along with element abundances w ere interpreted by 
IZubko. Krelowski fc Wegnerl l| 19961 .1 19981 ) using rather wide 
size distributions of grains. These authors applied spherical 
particles of different materials as well as inhomogeneous par- 
ticles, which makes impossible a direct comparison of their 
results with those obtained in our work. 

HD 197770. For this star the polarizing efficiency is max- 
imal: it reaches P/A = 3.7 at the J band (A = 1.25 fim). 
The polarization is enhanced in the region of the UV bump 
(super-Serkowski behaviour) but strong UV bump in extinc- 
tion results in no structure for P/A curve (see insert in 
Fig. |3} . Our model shows the largest fraction of graphite in 
the grain mixture for this star. 

Thus, we have fitted the observational data using pro- 
late or oblate spheroids with a/b — 2 or 3 and intermedi- 
ate values of the alignment parameter So = 0.3 — 0.5 fim. 
Contributions of carbonaceous and silicate grains to extinc- 
tion and polarization are nearly equal excluding the case of 
HD 62542 for which extinction is very peculiar. A contribu- 
tion of small graphite particles to extinction usually is less 
than 10%. The parameters of the size distributions for sili- 
cate and amorphous carbon grains differ from the standard 




Figure 4. Comparison of observed and modelled normalized extinction (left panel) and polarizing efficiency (right panel) for five stars. 
The solid curves show results of calculations. The model parameters are given in Table [4] 
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Table 4. Modelling results. 



Star 




amorphous 


carbon 


q 


/C Si 


astrosil 


9 


graphite 


shape; 
a/b 


<5o(/im) 


Q (deg) 


HD 24263 


0.48 


0.08 


0.40 


3.5 


0.47 


0.04 


0.15 


2.1 


0.05 


pro; 2 


0.3 


60 


HD 62542 


0.85 


0.07 


0.10 


2.7 


0.10 


0.005 


0.25 


2.7 


0.05 


pro; 2 


0.5 


10 


HD 99264 


0.44 


0.09 


0.30 


6.0 


0.50 


0.04 


0.18 


2.0 


0.06 


pro; 3 


0.5 


50 


HD 147165 


0.49 


0.06 


0.30 


2.2 


0.44 


0.04 


0.20 


2.2 


0.07 


pro; 2 


0.5 


54 


HD 147933 


0.50 


0.10 


0.25 


2.0 


0.47 


0.07 


0.20 


2.0 


0.03 


obi; 3 


0.3 


33 


HD 179406 


0.45 


0.10 


0.40 


6.0 


0.47 


0.04 


0.18 


2.5 


0.08 


pro; 3 


0.5 


19 


HD 197770 


0.43 


0.05 


0.20 


2.8 


0.46 


0.01 


0.17 


1.3 


0.11 


pro; 3 


0.5 


42 



»V,mm and ry, max are in ^tm 



MRN mixture (see Sect. 12. 2) . As usual, very small particles 
with rv,mm 0.01 fim are absent in the grain ensembles. Be- 
sides this, in many cases the slope of the n{rv) dependence 
is less steep than that for the MRN mixture (i.e. q < 3.5). 

The inclination of the magnetic field direction relative 
to the line of sight determined for all target stars lies in 
the interval Q ss 30 - 60°. These values of fi are derived 
rather accurately because of a strong dependence of P/A on 
Q (see Figs. 4 and 6 in VD08). The magnetic fields seem 
to be connected with local dust clouds in the directions to 
the selected stars. At the same time, the obtained values 
of f2 do not contradict the general pattern of the galactic 
magnetic field in the solar neighbourho od as found from 
pulsar measurements (|Valleell2004l . l2008l ). 

It should be noted that one usually assumed that small 
grains were badly oriented as the interstellar polariz ation de- 
gree decreases with a growing A -1 in the UV region |Whittetl 
l200St ). However, for all the objects we have fitted the ob- 
servational data by using the simple size and orientation 
distributions with small particles being rather well aligned. 
To see that one should consider Eqs. (|lip - (|13[l and keep in 
mind that <5o is about 0.4 and the temperature ratio about 
0.1. For ry ~ 0.06/im, this gives £ ~ 0.5 when £ = cor- 
responds to the perfect alignment while £ = 1 means the 
random orientation. 



4 CONCLUSIONS 

We have extended the spheroidal model of interstellar dust 
by considering imperfectly dynamically aligned amorphous 
carbon and silicate prolate and oblate spheroids with a 
power-law size distribution and a small fraction of graphite 
spheres involved to explain the 2175A bump. Our model be- 
ing as simple as possible was applied to fit the the normalized 
extinction A(X)/Av and the polarizing efficiency P(X)/A(X) 
observed in the near IR to far UV region. 

In contrast to earlier works we exactly calculated the 
optical properties of spheroids, for the first time utilized 
their imperfect Davis-Greenstein alignment and applied the 
spheroidal model to fit the data obtained for individual stars 
instead of the commonly used average Galactic curves. 

A comparison of our model with those recently sug- 
gested by other authors has demonstrated that the extinc- 
tion data are well fitted by all the models while only our 
model allows one to satisfactorily fit the polarization data 



with simple size and orientation distributions. Contrary to 
other works we find that the data can be well fitted with 
aligned carbonaceous and small silicate non-spherical grains. 

Investigation of our model has shown that the only pa- 
rameter that can be determined more or less accurately from 
observational data fitting is the direction of magnetic field 
relative to the line of sight. As the projection of the mag- 
netic field on the sky plane is given by the positional angle of 
polarization #, multi-wavelength extinction and polarization 
observations can be used for study of the spatial structure 
of interstellar magnetic fields. 

Variations of the interstellar extinction and polariza- 
tion curves for the considered stars can be interpreted by 
differences in the relative fractions of carbon and silicate 
grains and their size distributions in the clouds crossed by 
the lines of sight with the alignment parameter having in- 
termediate values So — 0.3 — 0.5/im (for the particle aspect 
ratio a/b — 2 — 3) . We assume that ambiguous results could 
be ruled out by considering dust phase abundances. 
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